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GENE RATING A NONLINEAR MODEL AND GENERATING 
DRIVE SIGNALS FOR SIMULATION TESTING USING 

THE SAME 

BACKGROUND OF THE INVENTION 

The present invention relates to computer 
modeling of a system, machine or process. More 
particularly, the present invention relates to 
generating an improved nonlinear system model, such as 
a neural network, as well as using a nonlinear system 
model or the improved nonlinear system model to 
generate drive signals as input to- a vibration system. 

Vibration systems that are capable of 
simulating loads and/or motions applied to test 
specimens are generally known. Vibration systems are 
widely used for performance evaluation, .durability 
tests, and various other purposes as they are highly 
effective in the development of products. For 
instance, it is quite common in the development of 
automobiles, motorcycles, or the like, to subject the 
vehicle or a substructure thereof to a laboratory 
environment that simulates operating conditions such 
as a road or test track. Physical simulation in the 
laboratory involves a well-known method of data 
acquisition and analysis in order to develop drive 
signals that can be applied to the vibration system to 
reproduce the operating environment. This method 
includes instrumenting the vehicle with transducers 
"remote" to the physical inputs of the operating 
environment. Common remote transducers include, but 
are not limited to, strain gauges, accelerometers, and 
displacement sensors, which implicitly define the 
operating environment of interest . The vehicle is then 
driven in the same operating environment, while remote 
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transducer responses (internal loads and/or motions) 
are -recorded. During simulation with the vehicle 
mounted to the vibration system, actuators of the 
vibration system are driven so as to reproduce the 
5 recorded remote transducer responses on the vehicle in 
the laboratory. 

However, before simulated testing can occur, 
the relationship between the input drive signals to 
the vibration system and the responses of the remote 

10 transducers must be characterized or modeled in the 
laboratory. This procedure is referred to as "system 
identification" and involves obtaining a respective 
linear model or transfer function of the complete 
physical system (e.g. vibration system, test specimen, 

15 and remote transducers) hereinafter referred to as the 
"physical system" and calculating an inverse system 
linear model or transfer function of the same. The 
inverse linear system model or transfer function is 
then used iteratively to obtain suitable drive signals 

20 for the vibration system to obtain substantially the 
same response from the remote transducers on the test 
specimen in the laboratory situation as was found in 
the operating environment. 

As those skilled in the art would 

25 appreciate, this process of obtaining suitable drive 
signals is not altered when the remote transducers are 
not physically remote from the test system inputs 
(e.g. the case where "remote" transducers are the 
feedback variables, such as force or motion, of the 

30 vibration system controller) . 

Although the above -described method for 
modeling and obtaining drive signals for a vibration 
system has enjoyed substantial success, there is a 
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continuing need to improve such systems. In 
particular, there is a need to improve the system 
model and/or the process for obtaining the drive 
signals . Current techniques for obtaining the drive 
5 signals subject the test specimen to loads and 
displacements during the iterative process to 
ascertain the drive signals. These initial loads and 
displacements can damage the test specimen before the 
intended testing occurs. Accordingly, a method and 
10 system that reduces application of "non- testing" loads 
and displacements in developing drive signals would be 
desirable. 

SUMMARY OF THE INVENTION 
A method and system for generating an 

15 improved nonlinear system model includes generating a 
linear system model and using a response therefrom to 
generate the nonlinear system model. In one 
embodiment, the method includes applying an original 
system input to the system and" ^obtaining a first 

2 0 system output from the system; generating a linear 
system model; applying an alternative input to the 
linear system model and obtaining a linear system 
output from the linear system model; applying the 
alternative input to the system and obtaining a second 

25 system output, and generating a nonlinear system model 
using the alternative input, the linear system output 
from the linear system model and the second output 
from the system. 

A method and system for generating drive 

30 signals for a test system uses the improved nonlinear 
system model or a conventional nonlinear system model. 
In a first exemplary embodiment, the method includes 
applying a drive signal to the system; obtaining a 
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response of the system to the drive signal; generating 
an inverse nonlinear system model using the drive 
signal and the response; and calculating a final drive 
signal using a selected response and the inverse 
5 nonlinear system model - 

In a second embodiment for generating drive 
signals, the method includes applying a" drive signal 
to the system; obtaining a response of the system to 
the drive signal; generating a nonlinear system model 
10 using the drive signal and the response; generating an 
inverse linear system model; calculating a test drive 
signal using a selected response and the inverse 
linear system model; applying the test drive signal to 
the nonlinear system model; obtaining an actual 
25 is response of the nonlinear system model to the test 

drive signal; calculating an error as a function of 
the actual response and the selected response; and if 
the error exceeds a selected threshold, obtaining a 
new drive signal using the linear system model and 
repeating the above -de scribed steps where the test 
drive signal comprises the new drive signal until the 
error is less than or equal to the selected threshold. 



30 

20 



BRIEF DESCRIPTION OF THE DRAWINGS 
25 FIG. 1 is a block diagram of an exemplary 

environment for practicing the present invention. 

FIG. 2 is a computer for implementing the 
present invention. 

FIG. 3A is a flow chart illustrating the 
30 steps involved in a system identification phase of a 
prior art method of vibration testing. 

FIG. 3B is a flow chart illustrating the 
steps involved in an iterative phase of a prior art 
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method of vibration testing. 

FIG. 3C is a flow chart illustrating the 
steps involved in another iterative phase of a prior 
art method of vibration testing. 
■ 5 FIG. 3D is a flow chart illustrating the 

steps involved in another iterative phase of vibration 
testing. 

FIG. 4A is a detailed block diagram of a 
prior art iterative process for obtaining drive 
10 signals for a vibration system. 

FIG. 4B is a detailed block diagram of 
another prior art iterative process for obtaining 
drive signals for a vibration system. 

FIG. 4C is a detailed block diagram of 
15 another iterative process for obtaining drive signals 
for the vibration system. 

FIG. 5 is a block diagram/flow chart of the 
prior art system identification/iterative process for 
obtaining drive signals for a vibrat-ion system. 

30 

20 FIG. 6 is a block diagram/flow chart of a 

system identification/iterative process for obtaining 
drive signals for a vibration system of the present 
invention. 

FIG. 7 is a flow chart illustrating the 
25 steps involved in a system identification phase of an 
exemplary method for obtaining a nonlinear system 
model . 

FIG. 8A is a detailed block diagram of a 
first iterative process for obtaining drive signals 
30 for a vibration system according to the present 
invention. 

FIG. 8B is a detailed block diagram of a 
second iterative process for obtaining drive signals 
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for a vibration system of the present invention. 

FIG. 8C is a detailed block diagram of a 

third iterative process for obtaining drive signals 

for a vibration system of the present invention. 
5 FIG. 9 is a block diagram/flow chart of a 

fourth process for obtaining drive signals for a 

vibration system according to the present invention. 

FIG. 10 is a portion of a flow chart 

illustrating a method in accordance with FIG. 9. 
0 FIG. 11 is a block diagram of training a 

neural network structure for an inverse nonlinear 

model . ~~ 

FIG. 12 is a block diagram of training a 

neural network structure for a forward nonlinear 
5 model . 

FIG. 13 is a flow chart illustrating a 
method for generating an improved nonlinear system 
model . 

FIG. 14 is a flow chart illustrating a 
0 second method for generating an improved nonlinear 
system model . 

DETAILED DESCRIPTION OF THE ILLUSTRATIVE EMBODIMENTS 
FIG. 1 illustrates a physical system 10. The 

5 physical system 10 generally includes a vibration 
system 13 comprising a servo controller 14 and an 
actuator 15. In the schematic illustration of FIG. 1, 
the actuator 15 represents one or more actuators that 
are coupled through a suitable mechanical interface 16 

0 to a test specimen 18. The servo controller 14 
provides an actuator command signal 19 to the actuator 
15, which in turn, excites the test specimen. 18. 
Suitable feedback 15A is provided from the actuator 15 
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to the servo controller 14. One or more remote 
transducers 20 on the test specimen 18, such as 
displacement sensors, strain gauges, accelerometers, 
or the like, provide a measured or actual response 21. 
5 A physical system controller 23 receives the actual 
response 21 as feedback to compute a drive 17 as input 
to the' physical system 10. In an iterative process 
discussed below, the physical system controller 23 
generates the drive 17* for the physical system 10 

10 based on the comparison of a desired response provided 
at 22 and the actual response 21 of the remote 
transducer 20 on the test specimen 18. Although 
illustrated in FIG. 1 for the single channel case, 
multiple channel embodiments with response 21 

15 comprising N response components and the drive 17 
comprising M drive components are typical and 
considered another embodiment of the present 
invention. 

Although described hereinwhere the physical 
20 system comprises the vibration system 13 and remote 
transducer 20, aspects of the present invention 
described below can be applied to other physical 
systems. For instance, in a manufacturing process, the 
physical system includes the manufacturing machines 
25 (e.g. presses, molding apparatus, forming machines, 
etc.) and the drive 17 provides command signals to 
said machines, and the actual response 21 comprises 
manual or automatic measured parameters of the 
manufactured article such as a critical dimension. 
3 0 Another example includes an oil refinery where the 
physical system is the process plant and the actual 
response 21 comprises intermediate or final parameters 
related to output products. 
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FIG. 2 and the related discussion provide a 
brief, general description of a suitable computing 
10 environment in which the invention may be implemented. 

Although not required, the physical system controller 
5 23 will be described, at least in part, in the general 
context of computer- executable instructions, such as 
15 program modules, being executed by a computer 30. 

Generally, program modules include routine programs, 
objects, components, data structures, etc., which 
10 perform particular tasks or implement particular 
abstract data types. The program modules are 

20 

illustrated below using block diagrams and flowcharts. 
Those skilled in the art can implement the block 
diagrams and flowcharts to computer- executable 
15 instructions. Moreover, those skilled in the art will 

25 

appreciate that the invention may be practiced with 
other computer system configurations, including multi- 
processor systems, networked personal computers, mini 
computers, main frame computers,' and the like. The 

30 

20 invention may also be practiced in distributed 
computing environments where tasks are performed by 
remote processing devices that are linked through a 
communications network. In a distributed computer 
35 environment, program modules may be located in both 

25 local and remote memory storage devices. 

The computer 30 illustrated in FIG. 2 
comprises a conventional personal or desktop computer 
40 having a central processing unit (CPU) 32, memory 34 

and a system bus 36, which couples various system 
30 components, including the memory 34 to the CPU 32. The 
system bus 36 may be any of several types of bus 
^ 5 structures including a memory bus or a memory 

controller, a peripheral bus, and a local bus using 
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any of a variety of bus architectures. The memory 34 
includes read only memory (ROM) and random access 
memory (RAM) . A basic input /output (BIOS) containing 
the basic routine that helps to transfer information 
5 between elements within the computer 30, such as 
during start-up, is stored in ROM. Storage devices 38, 
such as a hard disk, a floppy disk drive, an optical 
disk drive, etc., are coupled to the system bus 36 and 
are used for storage of programs and data. It should 
10 be appreciated by those skilled in the art that other 
types of computer readable media that are accessible 
by a computer, such as magnetic cassettes, flash 
memory cards, digital video disks, random access 
memories, read only memories, and the like, may also 
15 be used as storage devices. Commonly, programs are 
loaded into memory 34 from at least one of : the storage 
devices 38 with or without accompanying data. 

An input device 40 such as a keyboard, 
pointing device (mouse), or the like, allows the user 
20 to provide commands to the computer 30. A monitor 42 
or other type of output device is further connected to 
the system bus 36 via a suitable interface and 
provides feedback to the user. The desired response 22 
can be provided as an input to the computer 3 0 through 
25 a communications link, such as a modem, or through the 
removable media of the storage devices 38. The drive 
signals 17 are provided to the physical system 10 of 
FIG. 1 based on program modules executed by the 
computer 30 and through a suitable interface 44 
30 coupling the computer 30 to the vibration system 13. 
The interface 44 also receives the actual response 21. 

Before describing the present invention, it 
may also be helpful to review, in detail, a known 
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method for modeling the physical system 10 and 
obtaining the drive 17 to be applied thereto. Although 
described below with respect to a test vehicle, it 
should be understood that this prior art method and 
5 the present invention discussed below are not confined 
to testing only vehicles, but can be used on other 
types of test specimens and substructures or 
components thereof. In addition, the description is 
done assuming spectral analysis based modeling 

10 estimation and implementation though operations can be 
carried by several other linear mathematical 
techniques (e.g. Adaptive Filter- Tuning or Adaptive 
Inverse Control (AIC) type models, parametric 
regression techniques such as Auto Regressive 

15 Exogenous (ARX) and State Space types of models, or 
combinations thereof) . 

Referring to FIG. 3A, at step 52, the test 
vehicle is instrumented with the remote transducers 
.20. At step 54, the vehicle is subjected to the field 

20 operating environment of interest and the remote 
transducer responses are measured and recorded. For 
instance, the vehicle can be driven on a road or test 
track. The measured remote transducer responses, 
typically analog, are stored in the computer 30 in a 

25 digital format through analog- to-digital converters, 
as is commonly known. 

Next, in a system identification phase, the 
input /output linear system model of the physical 
system 10 is determined. With the vehicle mounted on 

30 the laboratory test rig, this procedure includes 
providing drive 17 as an input to the physical system 
10 and measuring the remote transducer response 21 as 
an output at step 56. The drive 17 used for linear 
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model estimation can be random "white noise" having 
frequency components over a selected bandwidth. At 
10 step 58, an estimate of the linear model of the 

physical system 10 is calculated based on the input 
5 drive applied and the remote transducer response 
obtained at step 56. In one embodiment, this is 
15 commonly known as the "frequency response function" 

(FRF) . Mathematically, the FRF is a N x M matrix 
wherein each element is a frequency dependent complex 
variable (gain and phase versus frequency) . The 
columns of the matrix correspond to the inputs, while 
the rows correspond to the outputs. As appreciated by 
those skilled in the art, the FRF may also be obtained 
directly from prior tests using the physical system 10 
25 is or other systems substantially similar to the physical 

system 10 

An inverse linear model H(f)" 1 is needed to 
determine the physical drive 17 as a function of the 
remote responses at step 60. As appreciated by those 
skilled in the art, the inverse linear model can be 
calculated directly. Also, the term "inverse" model as 
used herein includes a M x N "pseudo-inverse" model 
for a non-square N x M system. 

At this point in the prior art, the method 
25 enters an iterative phase, illustrated in FIGS. 3B and 
4A, to obtain drive 17 which produces actual response 
21 that ideally replicates the desired remote 
40 transducer response 22 (hereinafter "desired 

response"). The inverse physical system linear model 
H(f)' 1 is represented at 72, while physical system 
{vibration system, test vehicle, remote transducers 
and instrumentation) is represented at 10. Referring 
to FIG. 3B, at step 78, the inverse linear model 72 is 
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applied to a target response correction 77 in order to 
determine an initial drive 17 x x (t) . The target 
response correction 77 can be the desired response 22 
for the initial drive, though most often it is reduced 
5 by a relaxation gain factor 95 . The calculated drive 
17 x ; (t) from the inverse linear model 72 is then 
applied to the physical system 10 at step 80. The 
actual remote transducer response 21 (hereinafter 
"actual response") y x (t)'of the physical system 10 to 
10 the applied drive 17 x^t) is then obtained at step 86. 
If the complete physical system 10 is linear (allowing 
a relaxation gain 95 of unity) , then the initial drive 
17 Xj(t) could be used as the required drive. However, 
since physical systems are typically nonlinear, the 
15 correct drive 17 has to be arrived at by an iterative 
process. (As appreciated by those skilled in the art, 
drive 17 used in previous tests for a similar physical 
system may be used as the initial drive.) 

The iterative process involves recording the 
30 20 first actual response y^t) resulting from the initial 

drive x 1 (t) and comparing it with the desired response 
22 and calculating a response error 8 9 Ay x as the 
difference at step 88. (The first actual response 
signal y 1 (t) is provided at 87 in FIG. 4A.) The 

2 5 response error 8 9 Ay a is compared to a preselected 
threshold at step 90 and if the response error 89 
exceeds the threshold an iteration is performed. 
Specifically the response error 89 Ay r is reduced by 
the relaxation gain factor 95 to provide the new 
30 target response correction 77. In this embodiment, the 
inverse transfer function H(f) -1 is applied to the new 
target response correction 77 to create a drive 
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correction Ax 2 94 (step 91) that is added to the first 
drive x, (t) 17A to give a second drive Xj(t) 17 at step 
92. The iteration process (steps 80-92) is repeated 
until the response error 89 is brought down below the 
5 preselected threshold on all channels of the response. 
The last drive 17, which produced a response 21, that 
was within the predetermined threshold of the desired 
response 22, can then be used to perform specimen 
testing. 

0 As described, the response error 89 Ay is 

commonly reduced by the relaxation gain factor (or 
iteration gain) 95 to form the target response 
correction 77. The iteration gain 95 stabilizes the 
iterative process ■ and trades off rate -of -convergence 

5 against iteration overshoot. Furthermore, the 
iteration gain 95 minimizes the possibility that the 
test vehicle will be overloaded during the iteration 
process due to nonlinearities present in the physical 
system 10. As appreciated by those 'skilled in the art, 

0 an iteration gain can be applied to the drive 
correction 94 Ax and/or the response error 89 . It 
should be noted in FIG. 4A that storage devices 38 can 
be used to store the desired response 22, the actual 
responses 21 and previous drives 17A during the 

5 iterative process. Of course, memory 34 can also be 
used. Also, a dashed line 93 indicates that the 
inverse linear model 72 is an estimate of the inverse 
of the physical system 10. The block diagram of FIG. 
4A, as discussed above, can be implemented by those 

0 skilled in the art using commercially available 
software modules such as included with RPCIII™ from 
MTS Systems Corporation of Eden Prairie, Minnesota. 
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At this point, a modified method of the 
prior art for calculating the drive can also be 
discussed. The modified prior art method includes the 
steps of the identification phase illustrated in FIG. 
5 3A and many of the steps of the iterative phase 
illustrated in FIG. 3B. For convenience, the iterative 
steps of the modified method are illustrated in FIG . 
3C and the block diagram is illustrated in FIG. 4B. As 
illustrated in FIG. 43, the calculation of the target 
10 response correction 77 is identical. However, if the 
response error 8 9 between the actual response 21 and 
the desired response 22 is greater than a selected 
threshold, then the target response correction 77 is 
added to a previous target response 79A at step 97 to 
15 obtain a new target response 79 for the current 
iteration. The inverse linear model 72 iS; applied to 
the target response 79 to obtain the new drive 17. As 
illustrated in FIG. 4B, the iteration gain 95 can be 
used for the reasons discussed above-. 

30 

20 a further iterative method for calculating 

the drive is illustrated in FIGS. 3D and 4C. As 
illustrated in FIG. 4C, a linear drive 96 is obtained 
by applying the inverse linear model 72 to the desired 
response 22. The linear drive 96 is used as a 
25 reference for comparison similar to how the desired 
response 22 was used in the methods described above. 
Accordingly, the inverse linear model 72 is applied to 
the actual response 21 to obtain a drive estimate 98 
(step 101) . If a drive error 99 between the drive 
3 0 estimate 98 and the linear drive 96 is greater than a 
selected threshold (step 103) , then the drive 
correction 94 is added to a previous drive 17A at step 
105 to obtain a new drive 17 for the current 
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iteration. As illustrated in FIG. 4B, the iteration 
gain 95 can be used for the reasons discussed above. 
10 FIG. 5 is another block diagram flow chart 

illustrating the prior art technique of obtaining the 
5 drive signal for the system 10. In FIG. 5, like 
reference numbers have been used to identify similar 
elements as discussed above. w System identification" 
as described above with respect to FIG. 3A is shown 
generally at 100, while the iterative phase, described 
10 with respect to FIGS. 3B and 4A is illustrated 
generally at 102. The iterative phase 102 sequentially 
corrects the drive signal 17 and applies it to the 
system 10 until the response error 89 is within the 
predetermined threshold of the desired response 22. 
15 The final iterated drive signal is obtained at 104. Of 
course, a similar iterative phase corresponding to 
FIGS. 3C, 3D, 4B and 4C could also have been used to 
generate the final iterated drive signal 104 . 

Using the format of FIG. - 5, a first aspect 

30 

20 of the present invention is illustrated in FIG. 6. 
Generally, FIG. 6 illustrates a method 120 to control 
the system 10 to produce a selected response, and in 
particular, to generate a suitable drive signal to 
produce the selected response. The method 120 includes 
25 obtaining a nonlinear model 122 of the system 10 in a 
system identification 123 and using the nonlinear 
model 122, rather than the system 10, during an 
iterative phase 124 (executed within system controller 
23 of FIG. 1) to calculate or generate the final drive 
30 signal 104 intended to produce the desired response 
22. 

An exemplary method of system identification 
123 is illustrated in FIG. 7, which includes steps 52- 
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60 as discussed above wherein an inverse linear system 
model 72 (e.g-f inverse FRF) is obtained. System 
identification 123 however, also includes a step 126 
that comprises calculating or generating the nonlinear 
5 system model 122. In general, the nonlinear system 
model 122 comprises a tapped delay line or structure 
123 and at least one nonlinear regression structure 
125A. The tapped delay line 123 comprises a set of 
series connected delay elements 127 that provide a set 

10 of delayed inputs for each input drive signal 
(commonly digital values obtained at a selected 
sampling frequency). In FIG . 6, one input drive signal 
129 is illustrated although in many applications a 
plurality of individual drive signals are present, 

15 each drive signal being delayed by a corresponding 
tapped delay line to yield a corresponding set of 
delayed inputs. Typically, a nonlinear regression 
structure is provided for each output from the 
nonlinear model 122. In the embodiment illustrated, 

20 two regression structures 12 5A and 125B are shown and 
provide nonl inear mode 1 output s 1 3 1 A and 13 IB , 
respectively, but it should be understood that any 
number of nonlinear model outputs can be present. 
Likewise, although three delay elements 127 are shown 

25 any number of delay elements 127 can be used for each 
drive signal. 

In one embodiment, the nonlinear regression 
structure 125A and 125B each comprises a neural 
network. Regression structures 125A and 125B are means 

30 for regression modeling. In place of the neural 
network other forms of regression structures (or means 
for regression modeling) such as " nearest neighbors" , 
"support vector machines", "projection pursuit", 
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"multiple adaptive regression spines", "hinging 
hyperplanes" , "radial basis functions" "constrained 
topological mapping", etc. can also be used. 
Calculation of nonlinear regression models is well 
5 known in the art. For instance, software programs such 
as MATLAB, Neural Network Toolbox available from The 
Mathworks, Inc. of Natick, Massachusetts can be used 
to calculate a neural network based on the random 
drive and response of the system 10 to the drive 

10 obtained at step 56 in FIG. 7. 

Of course other drive signals besides a 
random drive can be applied to the system 10 used to 
calculate the inverse linear model 72 and/or the 
nonlinear system model 122. For instance, synthetic 

15 random signals having a prescribed power spectrum can 
be used. Random signals obtained from previous tests 
on the same or similar test system can also be used. 
Random signals obtained by use of preliminary system 
models (i.e. analytic/computer models or equations 

20 representing the dynamics of the test system) can also 
be used. Furthermore, multiple channel random input 
signals may be statistically correlated or 
uncorrelated (i.e. w orthogonal" ) . 

It should also be noted that the nonlinear 

25 model 122 need not correspond to the system 10 for 
which drive signals are ultimately intended. In other 
words, the nonlinear model 122 may be generated using 
drive and response signals from test system A, wherein 
the nonlinear model 122 may be then used to generate 

30 drive signals for test system B, which differs from 
test system A in specimen, test rig, or both. However, 
the greater the similarity between test systems A and 
B, the better the results, in general. 
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The nonlinear system model 122 is a w forward 
model" that models the system 10. The nonlinear system 
model 122 is then used during the iterative phase 124, 
rather than the system 10, to iterate to the drive 

5 signal which realizes a response within a 
predetermined threshold of the desired response 22 . 
FIG . 8A illustrates in block diagram form, the 
iterative phase 124 wherein a response error 89 and 
drive corrections 94 are used during the iterative 

0 process. If desired, the iterative phase 124 can be of 
the form where the response error 8 9 and target 
response corrections 77 are used ~ which is illustrated 
in FIG. 8B. Likewise, the iterative phase 124 can be 
of the form where the drive error 9 9 and drive 

5 corrections 94 are used, as illustrated in FIG. 8C. 

Referring back to FIG. 6, in many cases, the 
final drive 104 obtained at the end of the iterative 
phase 124 is sufficiently accurate to obtain 
substantially the desired respons'e/22 when applied to 

0 the system 10. However, if desired, conventional 
iteration can be performed as described above with 
respect to FIGS. 3B, 3C, 3D, 4A, 4B and 4C, which is 
generally illustrated at 130, to further refine the 
final drive signal 104. Likewise, if desired, the 

5 nonlinear system model 122 can be modified by further 
training, followed by more model iterations. 

As appreciated by those skilled in the art, 
the inventive method described above substantially 
reduces the number of drive signals applied to the 

0 system 10 in order to obtain the appropriate drive 
signal to generate the desired response 22 . 
Accordingly, wear and/or possible damage to the system 
10 is substantially reduced. 
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It should be noted that the inverse linear 
system model 72 can be obtained using any of the 
following techniques. As illustrated in FIGS. 6 and 7,- 
a forward linear model can be calculated using the 
5 same drive and response as those used for calculating 
the nonlinear system model 122, followed by inversion 
using conventional techniques. 

Alternatively, the inverse linear system 
model 72 can also be obtained by re- identification of 
10 the forward linear model, wherein a prescribed drive 
signal is applied to the nonlinear system model 122, 
obtaining the subsequent response, and then 
calculating the forward linear model by conventional 
techniques. The forward linear model can then be 
15 inverted using conventional techniques. This approach 
offers the possibility of using a drive .signal with 
characteristics different than those used for the 
nonlinear system model. 

In yet another method, direct calculation of 
20 the forward linear model can be obtained by 
linearization of the nonlinear model equations, 
followed by conventional inversion. 

The inverse linear model 72 can also be 
35 obtained directly without first calculating the 

25 forward linear model. For instance, calculation of the 
inverse linear model 72 can be performed directly 
using the same drive and response as those used for 
40 the nonlinear system model 122. In addition, the 

inverse linear model 72 can be re-identified by 
30 playing a prescribed drive signal into the nonlinear 
system model 122 and obtaining the subsequent 
45 response. The inverse linear model 72 can then be 

obtained by conventional techniques. Again, this 
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approach offers the possibility of using a drive 
signal with characteristics different than those used 
for obtaining the nonlinear system model 122, which 
may be desirable for some desired responses. Also, the 
5 inverse linear model 72 can be directly calculated by 
linearization of equations in an inverted nonlinear 
system model . 

The inverse linear model 72 described with 
respect to the method of FIG. 6 was static (i.e. non- 
10 changing values) . However, adjustments to the inverse 
linear model 72 can also be made during the iteration 
phase as described in co-pending application Serial 
No. 09/234,998, entitled, " METHOD AND APPARATUS FOR 
GENERATING DRIVE SIGNALS IN A VIBRATION TEST SYSTEM" , 
15 filed January 21, 1999, and incorporated herein by 
reference in its entirety. 

A second and separate aspect of the present 
invention is illustrated in FIGS 9 and 10 and also is 
a method 140 of controlling a test- system responsive 

30 

20 to a drive signal to produce a selected response. 
Referring to FIGS. 9 and 10, the method 140 includes a 
system identification process 144 (FIG. 9) wherein a 
test drive signal 142 (FIG. 9) is applied to the 
system 10 at step 150 (FIG. 10) . A subsequent response 
25 148 of the system 10 to the test drive signal 142 is 
then obtained at step 152. Using the test drive signal 
142 and subsequent response 148, an inverse nonlinear 
system model 154 (similar in structure to nonlinear 
model 122) is calculated at step 156. By applying the 
3 0 selected or desired response 22 to the inverse 
nonlinear system model 154, a final drive signal 158 
45 is calculated. In view that the inverse nonlinear 

system model 154 was calculated and used to obtain the 
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final drive signal 158, the model 154 more accurately 
models the system 10 than that which was obtainable 
from the linear system models, discussed above. In 
many cases, the final drive signal 158 is sufficiently 
accurate so that further iterations need not be 
necessary. However, if desired, iterative refinement 
of the drive signal 158 can be calculated using the 
inverse nonlinear system model 154 in a manner similar 
to that discussed above. 

It should also be noted that iterative 
processing can take the form of any of the techniques 
discussed above, with respect to -FIGS. 3B, 3C, 3D, 4A, 
4B and 4C. Specifically, drive corrections can be 
obtained from the inverse nonlinear system model or 
from the inverse linear model as used in the prior 
art. Alternatively, target response corrections can be 
applied to the inverse linear system model 154 wherein 
the new drive signal is obtained therefrom. 

In many instances, the ' ctesired response 22 
is prerecorded so the final drive 158 can be 
calculated by application of the inverse nonlinear 
system model 154 to the desired response 22 without 
regard to how long processing might take. The final 
drive 158 can then be stored and applied to the system 
10 as desired. However, in another embodiment, the 
nonlinear system model 154 can be suitably stored and 
accessed by a controller (e.g. computer or digital 
signal processor) for real time processing. In this 
manner, the desired response 22 can be an input to the 
controller wherein the controller calculates the drive 
for the system 10 using the nonlinear system model 154 
on a real time basis in order that the system 10 
subsequently produces the desired response as the 



WO 00/23934 



PCT/US99/24345 



-22- 

actual response. This allows testing of the test 
specimen without first recording a desired response. 

As indicated above, the nonlinear system 
model 122 and the inverse nonlinear system model 158 
can be generated using a number of nonlinear modeling 
techniques. One particularly useful technique is a 
neural network. 

FIG. 11 illustrates training of a neural 
network structure 160 for an inverse model. For 
purposes of simplification, a single input /single 
output embodiment is illustrated. A physical system 
response 162 is provided to a tapped delay line 164 
having a series of connected delay elements 166 that 
provide a plurality of delayed model inputs 168 to the 
neural network structure 160. A physical system drive 
172 is used as the model output. If desired, the 
output 172 can be provided to a tapped delay line 174 
having the series connected delay elements 166 to form 
a plurality of additional delayed model inputs 178 to 
the neural network structure 160. As appreciated by 
those skilled in the art, the number of delay elements 
166 present in the tapped delay line 164 and the 
tapped delay line 174 can be chosen independently. The 
neural network structure 160 is a well-known 
multilayer perceptron network type wherein the number 
of hidden layers, the number of neurons in each hidden 
layer, and the number of inputs and the number of 
outputs can be chosen independently. In addition, the 
form of the activation functions present in the neural 
network structure 160 can take well-known forms such 
as hyperbolic tangent and logistic sigmoid, to name a 
few. 

In the embodiment illustrated in FIG. 11, an 
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inverse neural network model is created wherein the 
model- input comprises a physical system response and 
the model output comprises a physical system drive. 
Forward neural network models such as used in the 

5 embodiment described above with respect to FIG. 6 can 
be similarly developed. FIG. 12 illustrates training 
of a forward neural network model wherein the input 
162 comprises the physical system drive and the output 
172 comprises the physical system response. 

0 As discussed above, multiple input/multiple 

output models can be created. Typically, a nonlinear 
system model, such as the neural network structure 
160, is formed for a multiple input/single output 
application wherein a complete neural network model 

5 having multiple outputs is formed from a plurality of 
neural network structures, each providing : a separate 
output as a function of a plurality of inputs. 

Although known methods for generating 
nonlinear models, such as generation of a neural 

0 network model as described above can be used in the 
aspects of the present invention discussed above, an 
improved method for nonlinear model generation is 
illustrated in FIG. 13. However, it should be 
understood that the method of FIG. 13 is not confined 

5 to use in modeling a vibration system in order to 
obtain suitable drive signals intended to produce" a 
desired response, but rather, can be used to generate 
a nonlinear model of any physical system (e.g. a 
manufacturing process, a complex machine, a component 

0 such as a shock absorber, etc.) comprising one or more 
inputs and one or more outputs. Accordingly, the 
method of FIG. 13 constitutes a third and separate 
broad aspect of the present invention. 
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In FIG. 13, a physical system 200 receives 
one or more original system input signals at 202 and 
1Q provides one or more first system output signals 204. 

The improved method of nonlinear model generation 
5 includes using the original system input signals 202 
and the first system output signals 204 in a 
conventional linear model generation process 205, such 
as a FRF model generation process to generate a linear 
system model 210 (in either a forward or inverse model 
10 application) . 

An alternative input signal (s) 209 is then 
applied to the linear system models 210 and a linear 
output prediction 212 (corresponding generally to the 
first output system signal (s) 204) is obtained 
15 therefrom. The alternative input signal (s) 209 is also 
applied to the physical system 200 and a second system 
output signal 211 is obtained. The alternative input 
signal (s) 209, the second system output signal 211 and 
the corresponding linear output' prediction 212 are 

30 

20 then used in a nonlinear model generation process 213 
(e.g. neural network modeling as discussed above) to 
generate a nonlinear system model 216 (in either a 
forward or inverse model application) . Unlike 
conventional nonlinear modeling generation, which 
25 forms a nonlinear system model as a function of only 
the alternative input signal (s) 209 and the second 
output signal 211, the nonlinear system model 216 
obtained is also a function of the corresponding 
linear response prediction 212. If the input signal (s) 
30 202 is subsequently applied to the nonlinear system 
model 216, a nonlinear response prediction 218 
accurately resembles the output signal 204, thereby 
verifying that the nonlinear system model 216 
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accurately models the physical system 200. 

In the embodiment illustrated, the separate 
input signals 202 and 209 have been used to generate 
the linear system model 210 and the inverse linear 
5 system model 216; however, as appreciated by those 
skilled in the art, the same input signals can be used 
This method is illustrated in FIG. 14. 

Although the present invention has been 
described with reference to preferred embodiments, 
10 workers skilled in the art will recognize that changes 
may be made in form and detail without departing from 
the spirit and scope of the invention. 
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WHAT IS CLAIMED IS: 

1. A method of controlling a test system 

responsive to a drive signal to produce a selected 
response , the method comprising : 

(a) applying a drive signal to the system; 

(b) obtaining a response of the system to 

the drive signal; 

(c) generating a nonlinear system model 
using the drive signal and the 
response ; 

(d> generating an inverse linear system 
model ; "~ 

(e) calculating a test drive signal using 

the selected response and the inverse 
linear system model; 

( f ) applying the test drive signal to the 

nonlinear system model; 

(g) obtaining an actual response of the 
nonlinear system model to the test 
drive signal; 

(h) calculating an error as a function of 

the actual response and the selected 

response; and 
if the error exceeds a selected threshold: 

(i) obtaining a new drive signal ; and 

(j) repeating steps (f) through (i) 
wherein the test drive signal is 
the new drive signal until the 
error is less than or equal to the 
selected threshold. 



2 . The method of claim 1 wherein obtaining a 

new drive comprises : 
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calculating a target response correction; 
calculating a drive correction using the 

inverse linear system model; and 
combining the drive correction with a 

previous test drive signal to obtain a 

new test drive signal. 



3 . The method of claim 1 wherein obtaining a 
new drive comprises: 

calculating a target response correction ; 
2Q combining the target response correction 

with a previous target response; and 
using the inverse linear system model to 

obtain the new drive. 

25 

4 . The method of claim 1 wherein generating the 
inverse linear system model comprises: 

generating a forward linear system model 
using the drive signal of step (a) and 

30 

the response signal of step (b) ; and 
calculating the inverse linear system model 
from the forward linear system model. 

35 

5. The method of claim 1 wherein generating the 
inverse linear system model comprises: 

applying a selected drive signal to the 
nonlinear system model ; 
40 obtaining a resultant response signal from 

the nonlinear system model; 
generating a forward linear system model 
using the selected drive signal and the 
45 resultant response signal; and 

calculating the inverse linear system model 
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f rom the forward linear system model . 

6 . The method of claim 1 wherein generating the 
inverse linear system model comprises : 

linearizing equations of the nonlinear 
system model to obtain a forward linear 
system model; and 

calculating the inverse linear system model 
from the forward linear system model. 

7. The method of claim 1 wherein generating the 
inverse linear system model comprises : 

generating the inverse linear system model 
using the drive signal of step (a) and 
the response signal of step (b) . 

8. The method of claim 1 wherein generating the 
inverse linear system model comprises : 

applying a selected drive signal to the 

nonlinear system model; 
obtaining a resultant response signal from 

the nonlinear system model; and 
generating the inverse linear system model 

using the selected drive signal and the 

resultant response signal. 

9. The method of claim 1 wherein generating the 
inverse linear system model comprises: 

linearizing equations of the inverse 
nonlinear system model to obtain the 
inverse linear system model. 



10. The method of claim 1 and further 
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comprising : 

calculating a linear system model and. a 
response from the linear system model; 
and 

wherein generating the inverse nonlinear 
system includes using the response from 
the linear system model. 

11. The method of claim 10 wherein obtaining the 
response from the linear system model includes 
applying an alternative input signal different than 
the first -mentioned drive signal. 

12 . A method of controlling a test system 
responsive to a drive signal to produce a selected 
response, the method comprising: 

(a) applying a drive signal to the system ; 

(b) obtaining a response of the system to 

the drive signal; 

(c) generating an inverse nonlinear system 
model using the drive signal and the 
response; and 

(d) calculating a final drive signal using 

the selected response and the inverse 
nonlinear system model. 

13. The method of claim 12 wherein generating 
the inverse nonlinear system model includes: 

training a neural network using the drive 
signal and the response. 

14 . The method of claim 12 and further 
comprising : 
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calculating a linear system model and a 
response from the linear system model; 
and 

wherein generating the inverse nonlinear 
system includes using the response from 
the linear system model. 

15. The method of claim 14 wherein obtaining the 
response from the linear system model includes 
applying an alternative input signal different than 
the first -mentioned drive signal. 

16. The method of claim 15 wherein the selected 
drive signal corresponds to the first-mentioned drive 
signal. 

17. The method of claim 15 wherein the selected 
drive signal does not correspond to the first - 
mentioned drive signal, and wherein generating the 
inverse nonlinear system model includes using the 
selected drive. 

18. A method of obtaining a non- linear system 
model of a system, the method comprising: 

(a) applying an original input to the system 

and obtaining a first system output 
from the system ; 

(b) generating a linear system model using 

the original input and the first system 
output ; 

(c) applying an alternative input to the 
linear system model and obtaining a 
linear system output from the linear 
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system model; and 
<d) applying the alternative input to the 

system and obtaining a second system 

output from the system; 
(e) generating a nonlinear system model 

using the alternative input, the second 

system output from the system and the 

linear system output from the linear 

system model . 

19. The method of claim 18 wherein the original 
input and the alternative input are the same. 

20. The method of claim 18 wherein the nonlinear 
system model comprises a forward nonlinear system 
model and wherein the inputs comprise drives and the 
outputs comprise responses. 

21 . The method of claim 18 wherein the nonlinear 
system model comprises an inverse nonlinear system 
model and wherein the inputs comprise responses and 
the outputs comprise drives. 

22. A vibration system for testing a test specimen, 
the vibration system comprising: 

an actuator system couplable to the test 
specimen ; 

a servo controller coupled to operate the 
actuator system; 

a system controller providing a drive signal to 
the servo controller, the system controller 
having means for generating, storing and 
using a nonlinear system model to generate 
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drive signals that when applied to the servo 
controller cause a desired response from the 
10 test specimen. 

23. The vibration system of claim 22 wherein the 
nonlinear system model comprises a forward nonlinear 

l 5 system model. 

24. The vibration system of claim 23 wherein the 
means for generating, storing and using the forward 
nonlinear system model to generate the drive signal 
comprises means for iteratively applying drive signals 
to the forward nonlinear system model. 
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25. The vibration system of claim 22 wherein the 
nonlinear system model comprises an inverse nonlinear 
system model.* 

26. The vibration system of claim 22 wherein the 
nonlinear system model comprises a regression 
structure and a tapped delay line providing a set of 
inputs to the regression structure. 



27. A computer readable medium including 

instructions readable by a computer, which when 
implemented, cause the computer to control a physical 
system responsive to a drive signal to produce a 
selected response, the instructions performing steps 
comprising: 

(a) applying a drive signal to the system; 

(b) obtaining a response of the system to 

the drive signal; 
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(c) generating a nonlinear system model 
using the drive signal and the 

10 response ; 

(d) generating an inverse linear system 
model ; 

(e) calculating a test drive signal using 
75 the selected response and the inverse 

linear system model; 

(f) applying the test drive signal to the 
nonlinear system model; 

(g) obtaining an actual response of the 
nonlinear system "model to the test 
drive signal; 

(h) calculating an error as a function of 
the actual response and the selected 
response; and 

if the error exceeds a selected threshold: 
(i) obtaining a new drive signal; and 
( j ) repeat ing steps- ( f ) through { i ) 
wherein the test drive signal is 
the new drive s ignal unt i 1 the 
error is less than or equal to the 
selected threshold. 

35 

28 . The computer readable medium of claim 27 and 

further comprising instructions for: 

calculating a linear system model and a 

40 

response from the linear system model; 
and 

wherein generating the inverse nonlinear 
system includes using the response from 

45 

the linear system model. 
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29. The computer readable medium of claim 28 

wherein obtaining the response from the linear system 
10 model includes applying an alternative input signal 

different than the first-mentioned drive signal. 
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30. A computer readable medium including 
instructions readable by a computer,' which when 
implemented, cause the computer to control a physical 
system responsive to a drive signal to produce a 
selected response, the instructions performing steps 
comprising : 

(a) applying a drive signal to the system; 

(b) obtaining a response of the system to 

the drive signal; * 

(c) generating an inverse nonlinear system 
model using the drive signal and the 
response; and 

(d) calculating a final drive signal using 

the selected response and the inverse 
nonlinear system model . 

31. The computer readable medium of claim 3 0 and 
further comprising instructions for: 

calculating a linear system model and a 
response from the linear system model; 
and 

wherein generating the inverse nonlinear 
system includes using the response from 
the linear system model. 

32. The computer readable medium of claim 31 
wherein obtaining the response from the linear system 
model includes applying an alternative input signal 
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different than the first -mentioned drive signal. 

33 . A computer readable medium including 

instructions readable by a computer, which when 
implemented, cause the computer to obtain a non- linear 
system model of a system, the instructions performing 
steps comprising: 

(a) applying an original input to the system 

and obtaining a first system output 
from the system; 

(b) generating a linear system model using 

the original input~and the first system 
output ; 

(c) applying an alternative input to the 
linear system model and obtaining a 
linear system output from .the linear 
system model ; and 

(d) applying the alternative input to the 
system and obtaining a second system 
output from the system; 

(e) generating a nonlinear system model 

using the alternative input, the second 
system output from the system and the 
linear system output from the linear 
system model . 



34. The computer readable medium of claim 33 

wherein the original input and the alternative input 
are the same. 



35. The computer readable medium of claim 33 

wherein the nonlinear system model comprises a forward 
nonlinear system model and wherein the inputs comprise 
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drives and the outputs comprise responses. 

10 36. The computer readable medium of claim 33 

wherein the nonlinear system model comprises an 
inverse nonlinear system model and wherein the inputs 
comprise responses and the 

15 



20 



25 



30 



35 



40 



45 



55 



WO 00/23934 



PCT/US99724345 



//// 



\ 
t 



V 



j 



UJ 

o 
o 



C2 

y3iioyiNOO 

IAI31SAS 



1 



CD 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT7US99/24345 



2/2/ 




CD 

LL 





tn 


AG 


LU 

O 00 




> ^ 




DE 







SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



3/et 



INSTRUMENT 
VEHICLE 



52 



RECORD VEHICLE 
RESPONSE DURING 
TEST EVENT 



54 



APPLY RANDOM DRIVE 
TO SYSTEM AND 
OBTAIN RESPONSE 



56 



58 



CALCULATE 
LINEAR MODEL 
H 



CALCULATE INVERSE 
LINEAR MODEL 




FIG. 3B 
OR 

FIG. 3C 



FIG. 3A 

PRIOR ART 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCTAJS99/24345 



OBTAIN FIRST DRIVE 
(Xt) 



ADD DRIVE CORRECTION 
TO PREVIOUS DRIVE 
<X M > 



78 



APPLY DRIVE (Xj)TO 
PHYSICAL SYSTEM 



80 



OBTAIN ACTUAL 
RESPONSE (Y. ) 



CALCULATE DRIVE 
CORRECTION USING H 



-1 



86 



88 



CALCULATE RESPONSE 
ERROR BETWEEN ACTUAL 
RESPONSE AND DESIRED 
RESPONSE 



91 



90 



IS RESPONSE ERROR 
GREATER THAN A 
SELECTED THRESHOLD? 



YES 



FIG. 3B 

PRIOR ART 



NO 



BEGIN TESTING 
USING LAST 
DRIVE 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



S/iJ 

78 



OBTAIN FIRST DRIVE 
(X,) 



80 



APPLY DRIVE (Xj)TO 
PHYSICAL SYSTEM 



CALCULATE NEW DRIVE 
USING H -1 



99 



OBTAIN ACTUAL 
RESPONSE (Yj) 



86 



88 



CALCULATE RESPONSE 
ERROR BETWEEN ACTUAL 
RESPONSE AND DESIRED 
RESPONSE 



97 



ADD TARGET RESPONSE 

CORRECTION TO 
PREVIOUS TARGET (Y M ) 



IS RESPONSE ERROR 
GREATER THAN A 
SELECTED THRESHOLD? 



YES 



FIG. 3C 

PRIOR ART 




NO 



BEGIN TESTING 
USING LAST 
DRIVE 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



OBTAIN FIRST DRIVE 
(X!) 



ADD DRIVE CORRECTION 
TO PREVIOUS DRIVE 

(*i-1> 



78 



APPLY DRIVE (Xj)TO 
PHYSICAL SYSTEM 



80 



OBTAIN ACTUAL 
RESPONSE (Y.) 



86 



101 



CALCULATE DRIVE 
ESTIMATE USING H 



CALCULATE DRIVE ERROR 
BETWEEN DRIVE 
ESTIMATE AND LINEAR 
DRIVE 



103 



105 



IS DRIVE ERROR 
GREATER THAN A 
SELECTED THRESHOLD? 



/YES 



FIG. 3D 



NO 



BEGIN TESTING 
USING LAST 
DRIVE 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 \ 



PCT/US99/24345 




SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



a/zt 




SUBSTITUTE SHEET (RULE 26) 




SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



/o/ai 



JO 



TEST DRIVE 
SIGNAL 



100- 




INVERSE 

LINEAR h- 
MODEL 

CALC 



TEST RESPONSE 
w SIGNAL 



(12 



DESIRED RESPONSE 



22 



1Q2 



4-I 



72 



ERROR 



INITIAL 
DRIVE 




h 



72^ 



ERROR 



2ND 
LAST 

CORRECTION/THRIVE 



- ^RROR „ 



FINAL 
ITERATED 
T DRIVE 



-fO 



2ND 




vDRIVE fc 









-fO 




FIG. 5 

PRIOR ART 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



///if 



TEST 
DRIVE SIGNAL 















NONLINEAR 
MODEL 
CALC 



TEST /23 
►•RESPONSE SIGNAL / ,„ 



DESIRED RESPONSE 
124 
\ 




FIG. 6 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24J45 



123 



INSTRUMENT 
VEHICLE 



52 



RECORD VEHICLE 
RESPONSE DURING 
TEST EVENT 



54 



J 



APPLY RANDOM DRIVE 
TO SYSTEM AND 
OBTAIN RESPONSE 



56 



58 



CALCULATE 
LINEAR MODEL 
H 



CALCULATE INVERSE 

LINEAR MODEL 
u-1 



60 



CALCULATE INVERSE 
NONLINEAR MODEL 



126 



FIG. 8A. 
FIG. 8B 

OR 
FIG. 8C 
SUBSTITUTE SHEET (RULE 26) 



FIG. 7 



* • 



WO 00/23934 PCT/US99/24345 





SUBSTITUTE SHEET (RULE 26) 



4 % 



WO 00/23934 PCT/US99V24345 



/5/i) 



<2 - 




SUBSTITUTE SHEET {RULE 26) 



WO 00/23934 



PCT/US99/2434S 



142 



■10 



DRIVE SIGNAL- 



144 



X 




,148 

1 



INVERSE 
NONLINEAR 
MODEL 
CALC 



RESPONSE SIGNAL 



r 



fS4 



INVERSE 

NONLINEAR 

MODEL 



DESIRED RESPONSE 



INVERSE 

NONLINEAR 

MODEL 



/SB 




70 



ERROR 



INVERSE 
NONLINEAR 
MODEL 



CORRECTION 



'10 



FIG. 9 

SUBSTITUTE SHEET (RULE 26) 



4 * 



WO 00/23934 



PCT/US99/24345 



17/it 



140 



APPLY DRIVE TO 
SYSTEM 



150 



152 



OBTAIN RESPONSE FROM 
SYSTEM 



156 



CALCULATE INVERSE 
NONLINEAR MODEL 



FIG. 10 



SUBSTITUTE SHEET (RULE 28) 



« » 



WO 00/23934 



PCT7US99/2434S 



/8/e/ 



162 



164- 



/74 



168 



ir " \' 



z A h 



178 



\ 


V. 


V 2 


Y U 
KM K-N 


U K-2 U K-1 










MULTILAYER 










PERCEPTRON 










NEURAL 










NETWORK 



FIG. 11 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99/24345 



/?/// 



/ « J 74 

t& £fU r /66 r 16 6 ft6 6 \ ^ r /6 6 ^ 



168 



/60 



1 whoi 



r Z* 



m 





U K -2 


Uk-M %M y K-2 Y K-1 






no OO MULTILAYER 






^^g^^T PERCEPTRON 






^^^W^ NEURAL 






<5^J^>^ NETWORK 



3 



FIG. 12 



SUBSTITUTE SHEET (RULE 25) 



WO 00/23934 



PCT/US99/24345 



20/21 



r 202 



r 



209 



2O0 



PHYSICAL 
SYSTEM 



XL 



s-205 



LINEAR MODEL 
GENERATION 
PROCESS 



LINEAR SYSTEM 
MODEL 



r ZOO 



PHYSICAL 
SYSTEM 



-2ft 



NONLINEAR MODEL 
GENERATION 
PROCESS 



202^ 



\Jr$6 sffl 



NONLINEAR SYSTEM 
MODEL 



NONLINEAR 
RESPONSE 
PREDICTION 



FIG. 13 



SUBSTITUTE SHEET (RULE 26) 



WO 00/23934 



PCT/US99Q4345 



21/2/ 



202 



200 



PHYSICAL 
SYSTEM 



204- 



205 



LINEAR MODEL 
GENERATION 
PROCESS 




210 



LINEAR SYSTEM 
MODEL 



LINEAR RESPONSE 
PREDICTION 



NONLINEAR MODEL 
GENERATION 
PROCESS 



•216 



NONLINEAR SYSTEM 
MODEL 



r 



218 



w NONLINEAR RESPONSE 
PREDICTION 



FIG. 14 



SUBSTITUTE SHEET (RULE 26) 



INTERNATIONAL SEARCH REPORT 



International application No. 
PCTAJS99/24345 



A. CLASSIFICATION OF SUBJECT MATTER 
IPC(6) :C06G 7/4B 
US CL : 703/8 

According to International Patent CUssiftctlkH) (IPC) or to both national clsjiification and IPC 



FIELDS SEARCHED 



Minimum documentation searched (classification system followed by classification symbols) 
U.S. : 703/8, 1, 2. 6, 7, 21; 318/568.22. 561; 212/275; 360/78.14; 701/97 



Documentation searched other than minimum documentation to the extent that such documents ire included in the Gelds searched 



Electronic dati base consulted during the international search (name of data base and, where practicable, search terms used) 
IEEE. STN 



DOCUMENTS CONSIDERED TO BE RELEVANT 



Category* 



Citation of document, with indication, where appropriate, of the relevant passages 



Relevant to claim No. 



A 
A 
A 
A 

A, P 



US 4,885,692 A (KURIHARA et al.) 05 December 1989, col. 3-8 
US 5,598,076 A (NEUBAUER et al.) 28 January 1997, col. 6-7 
US 5,602,689 A (KADLEC et al.) 11 February 1997, col. 7-76 
US 5,677,609 A (KHAN et al.) 14 October 1997, col. 4-8 
US 5,908,122 A (ROB1NETT et al.) 01 June 1999, col. 9-25 



1-36 
1-36 
1-36 
1-36 
1-36 



| | Further documents are listed in the continuation of Box C. | } See patent family annex. 



Special categories of CTtod documenti: 

document defining the general ttite of the aft which b net coniidtred 
to be of partioubtf ntrTHNi 

earner document yubtahod oa or after tb* international filing date 

document which nay throw doubti oa priority deim(i) or which » 
eked L> astabLkfa the publicalioa date of toother citation or other 
apectal reeaon (u specified) 

document referring lo an oral dbcloaurc. wc, exhmiuon or other 
oocuaeotpublobad prior to tha intsmaoortaJ filing dale but Ut*r thaa 



later document publisbad after tha international filing d*U or priority 
data and not b oaaflict with the application but cited to understand 
the principle or theory underlying tha invention 



document of panic alar relevance; the claimed invention cannrt be 
coDaidered novel or oanool be considered to inrohc an tn»emi»e ucp 
when the document b taken a loot 

document of particular relevance; tha claim ad invention cermoi be 
considered to involve an inventive »tcp when the document b 
combined wkfa one or more other racb document!, lucb combination 
being, obrioiu to a peraon drilled in the art 

doeunent nenber of the same patent family 



Date of the actual completion of the international search 
07 MARCH 2000 



Date of mailing of the international search report 

06APK200Q 




:VIN TESKA 
Telephone No. (703) 305-9704 



Name and Dialling address of the ISA/US 
Commixsjooer of Palcnti and Truiemtrks 
Box PCT 

WubJiiglon, D.C 20231 
Facsimile No. (703) 305-3230 



Form PCT/ISA/210 (second sheet)* July 1992)« 



